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Abstract 



We study asymptotics of the recurrence coefficients of orthogonal polynomials associated 
to the generalized Jacobi weight, which is a weight function with a finite number of algebraic 
singularities on [—1,1]. The recurrence coefficients can be written in terms of the solution of 
the corresponding Riemann-Hilbert problem for orthogonal polynomials. Using the steepest 
■ descent method of Deift and Zhou, we analyze the Riemann-Hilbert problem, and obtain 

complete asymptotic expansions of the recurrence coefficients. We will determine explicitly 
the order 1/n terms in the expansions. A critical step in the analysis of the Riemann-Hilbert 
problem will be the local analysis around the algebraic singularities, for which we use Bessel 
functions of appropriate order. 



> . 

1 Introduction 

o 

■ We consider the generalized Jacobi weight 

(N ■ n ° 

w (x) = (1 -x) a (l + xfh(x) Yl \x -x v \ 2K , for x £ (-1,1), (1.1) 

where no is a fixed number, with 

— 1 < x\ < X2 < ■ ■ ■ < x no < 1, 2A„ > — 1, A„ / 0, a,/5>— 1, 

^ ■ and with h real analytic and strictly positive on [—1,1]. The points x±,... ,x no are called the 

algebraic singularities of the weight. Throughout the paper we use xq = —1 and x no+ i = 1, for 
notational convenience. All the moments of w exist so that we have a sequence of orthogonal 
polynomials. Denote the n-th degree orthonormal polynomial with respect to the generalized 
Jacobi weight by p n (z) = r y n z ' n + • • • , where 7„ > 0. These orthonormal polynomials satisfy a 
three term recurrence relation 

xp n (x) = a n+1 p n+ i(x) + b n p n {x) + a n p n -i(x), 

and we will investigate the asymptotic behavior of the recurrence coefficients a n and b n as 
n — ► oo. The generalized Jacobi weight has been studied before from other points of view in 



H [Eg, |g, H among others. 

For the pure Jacobi weig ht (1 - x) a {\ + x)P exact expressions are known for the associated 
recurrence coefficients a n and. b n ^ see ^, [l!|. The asymptotic behavior is given by 

a n = — + 0(l/n 2 ), b n = 0(l/n 2 ), as n — > oo. 
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In a previous paper with Kuijlaars, Mclaughlin and Van Assche [16], we considered the modified 
Jacobi weig ht (1 -x) a (l +x) l3 h(x). There, we were able to obtain complete asymptotic expan- 
sions of the associated recurrence coefficients in powers of 1/n. It turned out that, as for the 
pure Jacobi weight, the order 1/n terms in the expansions vanished. The asymptotic behavior of 
the recurrence coefficients of orthogonal polynomials associated to the generalized Jacobi weight 
dL|) has b een studied before by Golinskii filf| . He has proven that 

a n = - + 0(l/n), b n = 0(l/n), as n — > oo. (1-2) 

In this paper we will give stronger asymptotics. We will prove that the 0(1/ n) terms in ( |1.2| ) 
can be developed into complete asymptotic expansions in powers of 1/n. Here, in contrast with 
the (modified) Jacobi weight, the order 1/n terms in the expansions will not vanish and we will 
determine an explicit expression for them. 

Our approach is based on the characterization of orthogonal polynomials via a Riemann- 
Hilbert problem, due to Fokas, Its and Kitaev |ll[, and on an application of the steepest descent 
method for Riemann-Hilbert problems of Deift and Zhou ||. We have already applied this 
technique to the modified Jacobi weig ht pi, and in our case, the general scheme is the same. 



The main difference lies in the fact that we now have to do a local analysis around the algebraic 
singularities as well (not just only around the endpoints ±1), which will be done with (modified) 
Bessel functions of appropriate order. In the present paper, we will emphasize the construction 
of the local parametrix near the algebraic singularities, which is new. It will turn out that the 
order 1/n terms in the expansions of the recurrence coefficients come from this parametrix. The 
Riemann-Hilbert approach has been applied before to orthogonal polynomials, see ||, [j], [||, P, 



14, 15, 16]. Our result is the following. 



Theorem 1.1 The recurrence coefficients a n andb n of orthogonal polynomials associated to the 
generalized Jacobi weight ( |l.l[) have a complete asymptotic expansion of the form 



1 | ^ A k (n) ^ B k (n) 
On~- + >^ k~ ' 6 «~>^ fc— > L3 

k=l k=l 



as n — > oo. The coefficients Ak(n) and B k (n) are explicitly computable for every k, and the 
coefficients with the 1/n term in the expansions are given by 



Ai(n) 



- \ v a/1 — x^ cos (2n arccos x v — Q u ) , (1-4) 



B\{n) = — \ v \Jl — xl cos((2n + 1) arccos x v — Q u ), 



where 



(1.5) 



»0 



»() 



a. 



arccos x v 

k=u+l J V k=l / 



a/1 - x 2 u f 1 ]ogh(t) dt 



x u f 

k T-i 



(1.6) 



The integral in (1.6) is a Cauchy principal value integral. 
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This theorem shows that n{2a n — 1) and nb n are oscillatory and asymptotically behave like 
a superposition of 7io WcivG functions A. v cos{uj v n — 4> u ) with amplitudes A v = —X u y / 1 — x 2 ,, 
frequencies tu v = 2 arccos x u , and phase shifts (b v which are different for n(2a n — 1) and nb n . The 
amplitude A v depends on the location and the strength of the singularity x u , while the frequency 
u v depends only on the location of x u . The strengths of the other singularities has influence on 
the phase shift cp u . This discussion shows that the 0(l/n) behavior of the recurrence coefficients 
is intimately related to the behavior of our weight near the singularities. Note that if we have 
no singularities (i.e. Ai = • • • = A no = 0) all the amplitudes in the wave functions vanish. This 
implies that the order 1/n terms in the expansions of the recurrence coefficients vanish, which 
is in agreement with the case of the modified Jacobi weight . 

Remark 1.2 We have restricted ourselves to determine only the the order 1/n terms in the 
expansions of the recurrence coefficients. It is possible to determine the higher order terms in 
the same way if we work hard enough, but the calculations will be a mess. 

We will now compare our result with a conjecture of Magnus [|17]] about the asymptotic 
behavior of the recurrence coefficients of orthogonal polynomials associated to the weight 

f B(l-x) a (l + x) l3 (x 1 -x) 2X , for x e (-l,xi), 

w(x) = < (17) 
[ A(l-x) a (l + x) /3 (x-a;i) 2A , for x S (xi, 1). 

where A and B are positive constants, with —1 < x\ < 1, where a, (3 > —1 and where 2A > — 1. 



This weight allows a jump at xi, and is of the form ( Llj ) only if A = B. The conjecture is the 
following. 



Conjecture of Magnus [17]. The recurrence coefficients of orthogonal polynomials associated 



to the weight (|1.7|) satisfy 

cos (^2n arccos x\ — 4/ilog(4n sin arccos x\) — +o(l/n), (l.f: 



1 M 

2 ~ ~n 
2M 



b n = cos^(2n + 1) arccosxi — 4^ log(4n sin arccos xi) — + o(l/n), (1-9) 

asm oo. Here 

H = i- log(S/A), M = ^y^T^^l-xl (1.10) 

$ = (a + A)7r - (a + + 2A) arccos x\ - 2 arg T(X + in) - arg(A + in). (1-H) 

We want to show that, as a consequence of Theorem [O], the conjecture is true for the case 
A = B. To this end, we need to reformulate the conjecture for this case. If A = B we have by 



( [LIOD and (|j~Tll ) that fj, = 0, M = (\\\/2)y/l - xj, and 

[ (a + X)ir — (a + (3 + 2A) arccos x\ , if A > 0, 
$ = < 

I (a + X)ir — (a + (3 + 2A) arccos x\ — 3ir, if A < 0. 
Inserting this into ( |l.g| ) and ( |1.9[) the conjecture becomes. 

Conjecture of Magnus for the case A = B. The recurrence coefficients of orthogonal 
polynomials associated to the weig ht A(l - x) a (l + xf\x - xi| 2A satisfy 

1 A 



1 — x? cos( 2n arccos x\ — <5 I + o(l/n), 
2 2n v 1 ' 



b n = — — yjl — x\ cos^(2n + 1) arccos x\ — + o(l/n), 



A 

n 

as n — > co, where $ = (a + A)7r — (a + (3 + 2A) arccos xi. 
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If we apply Theorem 



1.1 



for the case uq = 1 and h = A, and using the fact that, see [O ], 



1 

i Vl -t 2 t - x x 



we see that the conjecture is true for the case A = B. We even have more since we were able to 
obtain complete asymptotic expansions of the recurrence coefficients, allow the analytic factor 
h, and allow more singularities. The full conjecture (i.e. A ^ B) remains open. 

We will also compare our result with Nevai's result flH] for an even positive weight p(x)\x\ 2X 



on [—1, 1] with p and p' continuous and with 2A > —1. Nevai showed in [18] that the recurrence 



coefficient a n of orthogonal polynomials associated to this weight satisfies, 

a n = ^ + (-l) n+1 -^ + o(l/n), asn^oo. (1.12) 
2 2n 



We apply Theorem 1.1 to the weight w(x) = (1 — x 2 ) a h(x)\x\ 2Xl , with h even. So, w is of the 



form of Nevai's weight. Using the fact that h is even and x\ = 0, we have 



log h(t) dt _ 



-i Vl - t 2 t - xx 

since the integrand is an odd function. Since no = 1, a = (3, and arccosxi = tt/2 this implies 



by Q1.6|) that the phase constant vanishes. Prom (|1,3|) and (1.4) we then have 



^ + (-l) n+1 ^ + 0(l/n 2 ), asrwoo. 
2 2n 



This is in agreement with Nevai's result, see (1.12). The error is stronger since we have 0(l/n 2 ) 
instead of o(l/n). However, here we are dealing with an even weight of the form p(x)\x\ 2X with 
p analytic, and Nevai's weights also include cases where p is non-analytic. Nevai also showed 
that the error is 0(l/n 2 ) if p is constant. Note, since the phase constant <&i vanishes, that by 
( |1.5| ) the order 1/n term in the expansion of b n vanishes. This is in agreement with the fact that 
b n = for an even weight. 

The present paper is organized as follows. In Section 2 we formulate the theory of orthogonal 
polynomials as a Riemann-Hilbert (RH) problem for Y. In Section 3 we do the asymptotic 
analysis of this RH problem. There, we want to obtain, via a series of transformations Y 
T i — > S i — > -R, a RH problem for R with a jump matrix close to the identity matrix. Then, R 
is also close to the identity matrix. For the last transformation S i— > R we have to do a local 
analysis near the endpoints and near the algebraic singularities. The local analysis near the 



endpoints has already been done in [16], but near the algebraic singularities it is new and will 
be done in Section 4. In the last section we determine a complete asymptotic expansion of the 
jump matrix for R. As a result, we obtain a complete asymptotic expansion of R, which will be 



used to prove Theorem 1.1 



2 RH problem for Y 

In this section we will characterize the orthogonal polynomials via a 2 x 2 matrix valued RH 
problem. This characterization is due to Fokas, Its and Kitaev fIT|] . We will also write down 
the recurrence coefficients a n and b n in terms of the solution of this RH problem. 

We seek a 2 x 2 matrix valued function Y(z) = Y(z; n, w) that satisfies the following RH 
problem. 
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RH problem for Y: 

(a) Y{z) is analytic for z G C \ [—1,1]. 

(b) Y possesses continuous boundary values for x G (—1, 1) \ {xi, . . . ,x no } denoted by Y + (x) 
and Y_(x), where Y+(x) and YL(x) denote the limiting values of Y(z') as z' approaches x 
from above and below, respectively, and 



y + {x) = y_(.x) 



1 w(x)\ r , . 1U , -. 

l J ' for x G (—1, 1) \ {xi, . . . , x no }. 



(c) Y{z) has the following asymptotic behavior at infinity: 

' 



Y(z) = (I + 0(l/z)) 



z~ 



as z — > oo. 



(d) Y{z) has the following behavior near z = 1: 

'1 \z-l\ a 



Y(z) 



O 



1 \z-l\ c 



if a < 0, 



< 5!::!! i. «°=* 



o 



i i 
i i 



if a > 0, 



as z -» 1, z G C\ [-1,1]. 
(e) l"(z) has the following behavior near z = — 1: 

'1 \z + l\P 



Y(z) 



O 



o 



1 \z + lf 

1 log|z+l| 
1 logjz+lj 

o' 1 1 



if p < 0, 
if /? = 0, 
if /? > 0, 



1 1 

as z -1, z G C\ [-1,1]. 
(f) 1" (z) has the following behavior near z = x u , for every v = 1, . . . , uq: 



Y(z) 



|2A t 



o 



i i 
i i) ' 



if K > 0, 



as z — ► £„, z G C \ [—1, 1] 

Remark 2.1 The O-terms in 
|z-l| 
|z-l| 



(2.1) 



(2.2) 



(2.3) 



(2.4) 



(2.5) 



( |2.4| ) and ( |2.5| ) are to be taken entrywise. So for example 
Y(z) = O I ! |* ~ |£ ] means that F u (z) = O(l), Yi 2 (z) = 0(|z - l| a ), etc. 



If we take care of the algebraic singularities x v of the generalized Jacobi weight in the same 
way as of the endpoints ±1 in [l6|, Section 2] we obtain the following theorem. 
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Theorem 2.2 The RH problem for Y has a unique solution Y{z) = Y(z;n,w) given by, 



y(z) 



(2.6) 



where ix n is the monic polynomial of degree n orthogonal with respect to the weight w and with 
7 n the leading coefficient of the orthonormal polynomial p n . 



The recurrence coefficients a n and b n can be written in terms of Y, see 
known || that 

a 2 = lim z 2 Yi 2 (z;n, w)Y 2 i(z;n,w), 



ni, m. it is 

(2.7) 



b n = lim (z — Yii(z;n + l,w)Y 2 2(z;n,w)). (2.8) 

So, in order to determine the asymptotics of the recurrence coefficients, we need to do an 
asymptotic analysis of the RH problem for Y. 



3 Asymptotic analysis of the RH problem for Y 

In this section we will do the asymptotic analysis of the RH problem for Y. The idea is to 
obtain, via a series of transformations 

Y^T^S^R, 

a RH problem for R which is normalized at infinity (i.e. R(z) — > / as z — > oo) and whose jump 
matrix is close to the identity matrix. As a result, the solution of the RH problem for R is also 
close to the identity matrix, cf. || ||. 

As mentioned in the introduction, we point out that the asymptotic analysis is analogous as 
in the case of the modified Jacobi weight, see Jl| . The main differences, which come from the 
algebraic singularities, are: 

• In every step we have to take care of the growth condition near the algebraic singularities, 
which was included in the RH problem for Y to control the behavior near these points. 

• In the second transformation T — > S the lens will be opened going through the algebraic 
singularities. 

• We have to do a local analysis around the algebraic singularities, not just only around the 
endpoints. This is a new and critical step in the analysis of the RH problem for Y, and is 
the most important difference with the case of the modified Jacobi weight. To emphasize 
this, the construction of the parametrix near the algebraic singularities will be done in a 
separate section. 



3.1 First transformation Y —>■ T 

We will first transform the RH problem for Y into a RH problem for T whose solution is bounded 
at infinity, and whose jump matrix has oscillatory diagonal entries. Let (f(z) = z + (z 2 — l) 1 / 2 
be the conformal mapping that maps C \ [— 1, 1] onto the exterior of the unit circle, and define 

T(z) = 2 n<T3 Y{z)y{z)- na ^ for z G C\ [-1,1], (3.1) 

where 03 = (J 2 %) is the Pauli matrix. Then, T is the unique solution of the following equivalent 
RH problem, cf. |L6], Section 3]. 
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RH problem for T: 

(a) T(z) is analytic for z G C \ [—1, 1]. 

(b) T(z) satisfies the following jump relation on (—1, 1) \ {x\, . . . , x„ }: 

T + (x) = T^(x)( ip+{x ^ 2n r ^l 2n ), ioTxe(-l,l)\{ Xl ,...,x no }. (3.2) 



P-(a?)~ 

(c) T{z) has the following behavior at infinity: 

T(z) = l + 0(l/z), asz^oo. (3.3) 

(d) T(z) has the same behavior as Y(z) as z — > 1, given by ( |2.3[ ), 

(e) T(z) has the same behavior as Y{z) as z — ► —1, given by ( |2.4| ). 

(f) T(z) has the same behavior as Y(z) as z — > av, given by (p.5|), for every i/ = 1, . . . , no- 



Remark 3.1 Condition (c) states that the RH problem for T is normalized at infinity Since 
l9 9 ±( a; )| = 1 f° r x G ( — 1) 1) we have by (|3.2| ) oscillatory diagonal entries in the jump matrix for 
T. 

3.2 Second transformation T — ► S 

We use the steepest descent method for RH problems of Deift and Zhou Q to remove the 
oscillatory behavior in ( |3.2[ ). See []|, ^] for an introduction. The idea is to deform the contour so 
that the oscillatory diagonal entries in the jump matrix for T are transformed into exponentially 
decaying off-diagonal entries. We then arrive at an equivalent RH problem for S on a lens shaped 
contour, with jumps that converge to the identity matrix on the lips of the lens, as n —* oo. 
This step is referred to as the opening of the lens. 

Since h is real analytic and strictly positive on [—1,1], there is a neighborhood U of [—1,1] 
so that h has an analytic continuation to U, and so that the real part of h is strictly positive 
on U. Hence, the factor (1 — x) a (l + x)^h(x) has a non- vanishing analytic continuation to 
z G *7\((-oo,-l] U [l,oo)), given by 

(l-z) a (l + zfh(z), 

with principal branches of powers. 

To continuate the factor \x— x„\ 2Xv analytically, where v G {1, . . . , no}, we divide the complex 
plane into two regions, which we denote by K l and , separated by a contour Y Xv going 
through x u , see Figure [j]. Here, K l Xv and K r Xv are the sets of all points on the left, respectively 
right, of Y Xv . We choose the contour T Xv so that the images of T Xi/ n C + and Y Xv n C_, 
under the mapping tp, are the straight rays, restricted to the exterior of the unit circle, with 
arguments arccosx*, and — arccosxjy, respectively. Here, C + is used to denote the upper half- 
plane {z | Imz > 0}, and C_ to denote the lower half-plane {z \ Imz < 0}. It turns out that T X[/ 
is a hyperbola and goes vertically through x v . We have made an exact plot of Y Xv for the case 
x v = 1/2, see Figure |l|. For v = 1, . . . , no, the factor \x — x u \ 2Xv has an analytic continuation to 
z G C \ Y Xv , given by 

J (x v -z) 2X ", for z G K Xv , 

\ {z-x u ) 2X % for z€KZ u , 

with again principal branches of powers. 
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1.6 - 








Figure 1: The contour T Xl/ for the case x u = 1/2. The star marks x u , and K x and K% are the 
sets of all points on the left, respectively right, of T Xu . 



Remark 3.2 It seems a bit awkward to work with this choice of r r „ instead of with the vertical 



line going through x u , but in Section |4.2| this will become clear. 



As a result, the generalized Jacobi weight w, given by (1.1), has a non- vanishing analytic 
continuation to z € U\ ((— oo, — 1] U [1, oo) U [U^^r^]), also denoted by w, given by 

w{z) = {I- z) a {l + zfh{z) 

v no 



x U u 

k=l 



where v = 0, . . . , no, and K, 



X'O 



K l 



n 

l=u+l 

c. 



\2Ai 



if z £ KZ. n K: 



(3.4) 



Remark 3.3 Note that only if \ u G N the analytic continuation of our weight is also analytic 
across the contour Y Xv . 



The jump matrix ( |3.2|) for T has the following factorization into a product of three matrices, 
based on the fact that ip + {x)if- (x) = 1 for x £ (—1, 1), 

/ ip + (x)~ 2n w(x) 







ip-(x)- 2n 

1 



w(x) 1 tp^(x 





-2n -j^ 



w(x) 





1 

w(xy 1 ip + (x)~ 2n i 



(3.5) 





-w(x)^ 1 

We note that 1/w does not have an analytic extension to a full neighborhood of (—1, 1). Instead, 
it has an analytic continuation to a neighborhood of (x u , x^+i), for every v = 0, . . . ,uq, where 
xq = —1 and x no _|_i = 1. We thus transform the RH problem for T into a RH problem for 
S with jumps on the oriented contour S, shown in Figure ||, that goes through the algebraic 
singularities x v . The precise form of the lens £ will be determined in Section 4.2. Of course it 
will be contained in U. We write 

S° = S \ { 1, x%i ■ ■ ■ , x no , 1}. 



Let us define, as in |16|, Section 4], 



S(z) 



T(z) 



T(z) 



-w(z) 1 ip(z) 

1 

w{z)- l tp{z)- 2n 



2n 



for z outside the lens, 

for z in the upper parts of the lens, 

for z in the lower parts of the lens. 



(3.6) 



Then, S is the unique solution of the following equivalent RH problem, cf. [16, Section 4]. 



S 



Figure 2: The contour S. 



RH problem for S: 

(a) S(z) is analytic for z G C \ S. 

(b) S'(z) satisfies the following jump relations on S°: 



S + (z) =< S_(z) 
5+(z) = S-(x) 



1 







™(z)-V0z)~ 2n i 





(c) 
(d) 



-w(x) 1 

S'(z) has the following behavior at infinity: 
S(z) = I + 0(1/ z), as z — > oo. 

For a < 0, S^ 2 ) nas the following behavior as z -> 1: 

S(*) = o(j |*ljja), asz^l,zG<C\£. 

For a = 0, «S(z) has the following behavior as z — ► 1: 

'fog | 2 — 1 1 log | Z — 1 1 



for z G S° n (C+ U C- 
for x G S°n (-1,1). 



S(z) = O 



log | Z — 1 1 log 1 2 — 1 1 

For a > 0, S'(z) has the following behavior as z — > 1: 



as z — ► 1, z G C \ S. 



S(z) 



O 






G 


r 






\z 


- 1| 


—a 


\z 


- ii 


—a 



as z — > 1 from outside the lens, 



, as z — > 1 from inside the lens. 



(3.7) 
(3.8) 

(3.9) 

(3.10) 

(3.11) 



(3.12) 



S(z) has the same behavior near —1 if we replace in (|3~10|) , (gH| ), and (|3l2|) , a by /?, 
|z — 1| by |z + 1| and take the limit z — » — 1 instead of z — » 1. 



(e) 

(f) For = 1, . . . , no, S(z) has the following behavior For A„ < we have 



S(z) = Ol 



1 |z — x u 



z — X, 



2K 
|2A„ I ' 



For X u > we have 



S(z) 



i r 
i in 



clS 2> ^ X i/ from outside the lens, 



(3.13) 



I £ I ~ 2A^ -j^ 

O ( , i-2A y ) > as z ~^ x v from inside the lens. 



(3.14) 
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Since |^(z)| > 1 for z £ C\ [—1, 1] we see from (|3.7f ) that the oscillatory terms on the diagonal 
entries in the jump matrix for T have been transformed into exponentially decaying off-diagonal 
entries in the jump matrix for S on the lips of the lens. So, the jump matrix for S converges 
exponentially fast to the identity matrix on the lips of the lens, as n — > oo. Hence, we expect 
that the leading order asymptotics are determined by the solution of the following RH problem. 

RH problem for N: 

(a) N(z) is analytic for z £ C \ [—1,1]. 

(b) N(z) satisfies the following jump relation on the interval (—1,1) \ {x±, . . . , x nQ }: 

w(x) s 



N+(x) =N-(x) [ w ^-i fovxe(-l,l)\{x 1 ,...,x no }. 

(c) N(z) has the following behavior at infinity: 
N(z) = l + 0(l/z), as z ^ oo. 



(3.15) 



(3.16) 



The solution of the RH problem for N is referred to as the parametrix for the outside region 
and it has been solved in [|l(], Section 5] using the Szego function associated with the generalized 
Jacobi weight w, 

n , ,_ {z-lYl 2 {z + iylZYKU(z-Xv) Xv / V-l) 1 / 2 f 1 iQgh(x) dx \ 

The Szego function D(z) associated to w is analytic and non-zero for z £ C\ [—1, 1], satisfies the 
jump condition D + {x)D-(x) = w(x) for x £ (—1, 1) \ {xi, . . . ,x no }, and -Dqo = lim 2 ^oo D(z) £ 



(0, +oo). The solution of the RH problem for N is then given by, see [16, Section 5], 



a(z)+a(z) 1 a(z)—a(z) 



N(z) = £>£ 



2i 



a(z)—a(z) 1 a(z)+a(z) 1 



D(zY 



-2i 



where 



a(z) 



(Z - 1)V4 



(z + iy/ 4 ' 

For later use we have the following lemma. 
Lemma 3.4 For every v = 0, . . . , no, 

D+(x) = y / w(x)e~ t ^ v( - x \ for x v < x < x v+ \. 
Here xq = —1, x no+ \ = 1, and 



no 



no 



a + arccos : 

v k=u+l J V k=l / 



y/l - x' 2 f 1 log h(t) dt 



-x 2 f 1 

k T-i 



VT^tft-x 

The integral in ( 3.21| ) is a Cauchy principal value integral. 



(3.18) 



(3.19) 



(3.20) 



(3.21) 
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Proof. We rewrite the expression ( [ml) for the Szego function as 



where 

1 f 1 logh(x) dt 



mz) 



2m y_! vT^t 5 * 



Now, we determine -D+(x) for x G (:r„, a^+i). So, we need to take the + boundary values for all 
quantities in ( 3.22| ). Using the Sokhotskii-Plemelj formula 12, Section 4.2] we have 



_ , . logh(x) 1 f 1 log h(t) dt 
$+(x) = — = + ' 



2tt« /_! 



2^1^^ 2ni J_ ly /T~fit 



where the integral is a Cauchy principal value integral, so that, by ( 3.22j ) and the fact that 



<p+(x) = exp(z arccos x), the lemma is proved after an easy calculation. □ 

Before we can to do the third transformation we have to be careful, since the jump matrices 
for S and N are not uniformly close to each other near the endpoints ±1 and near the algebraic 
singularities x v . Therefore, a local analysis near these points is necessary. Near the endpoints 



this has already been done in [16, Section 6]. 

We have constructed in [16, Section 6] a parametrix Pi in the disk Ugi with radius 5 > 0, 
sufficiently small, and center 1. This is a matrix valued function in U$ i, that has the same 
jumps as S on E, that matches with N on the boundary dUgi of Ugi, 

P 1 (z)N~ 1 (z) =I + 0(l/n), as n -> oo, uniformly for z G 8U s ,i, (3.23) 

and that has the same behavior as S(z) near z = 1. The parametrix Pi is given in 16, Section 



6], and is constructed out of Bessel function of order a. We note that the scalar function W in 
[fl6| , (6.27)], because of the extra factor Yiv=i \x — x u \ 2Xv in the generalized Jacobi weight, should 
have an extra factor n™li( z — x u ) 2Xv . 

Similarly we have constructed in [16, Section 6] a parametrix P_i in the disk Us-i with 
radius 5 > and center —1. This is a matrix valued function in Us-i that has the same jumps 
as 5 on E, that matches with TV on dUs -i 

P- 1 (z)N' 1 (z) =1 + 0(l/n), as n -» oo, uniformly for z G dU s ,-i, (3.24) 



and that has the same behavior as S(z) near z = —1. The parametrix is given in [16, Section 
6], and is constructed out of Bessel functions of order (3. We note that the scalar function W in 
PI, (6.52)] should have an extra factor n™Li(^ ~ z f K ■ 

We also have to construct a local parametrix P Xv near the algebraic singularities x v . Let Ug Xl/ 
be the disk, with center x v and radius 5 > so that the closures of the disks Us tX0 , . . . , Ug jXno+1 
don't intersect and so that all the disks lie in U. The construction of the parametrix P Xv will 
be done in Section |j. For now, let us assume that we have a 2 x 2 matrix valued function P Xv 
with the same jumps as S, that matches with N on dU$ jXv , 

P Xl ,(z)N- 1 (z)=I + 0(l/n), as n -> oo, uniformly for z G dU s , Xv , (3.25) 

and that has the same behavior as S(z) near z = x v . 
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Figure 3: The reduced system of contours with circles Ug Xv of radius 5 and center x v . 
3.3 Third transformation S — > i? 

Using the parametrix for the outside region, the parametrices near the endpoints, and the 
parametrices near the algebraic singularities we do the final transformation. Let us define the 
matrix valued function R as 

/ S(z)N^(z), for z 6 C \ (S U [U^+^J), 

= _ (3.26) 

[ 5(z)P 3 . i/ 1 (2;), for z € U s , Xu \ S, and 2/ = 0, ... ,n + 1. 

Remark 3.5 Note that the inverses of the parametrices exist. For N,P-\ and Pi this was 



already known, see [16]. In the next section we will show that P Xv is also invertible. 



If we take care of the behavior near the algebraic singularities in the same way as near the 
endpoints, it turns out that R satisfies the following RH problem, cf. |[(], Section 7], with jumps 
on the reduced system of contours £r, see Figure [|. 

RH problem for R: 

(a) R(z) is analytic for z £ C \ Y,r. 

(b) R(z) satisfies the following jump relations on E^: 

R + {z) = R^(z)P Xv (z)N~ 1 (z), 

for z € dUs iX „, and v = 0, . . . , no + 1, (3.27) 

1 0\ 



toxzeZ R \(u^ dUsxv y (3.28) 

(c) R{z) has the following behavior at infinity: 

R(z)=I + 0(l/z), asz^oo. (3.29) 



By (ITU and ( ^25|) , the jump matrices on the circles are uniformly close to the 



identity matrix as n — > cxd. On the lips of the lens we have by (3.28), as in [16, section 7], that 



the jump matrix converges uniformly to the identity matrix at an exponential rate. So, all jump 
matrices are uniformly close to the identity matrix. This implies that, cf. ||. @], 

R(z) = I + 0(l/n), as n — > 00, uniformly for z G C \ T,r. (3.30) 



Remark 3.6 We will show in Section 5 that the 0(l/n) term in ( 3.3C| ) can be developed into a 



complete asymptotic expansion in powers of 1/n. This expansion will be used to prove Theorem 



1.1 
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4 Parametrix near the algebraic singularity x v 

Fix i/G{l,...,rao}. In this section we construct a 2 x 2 matrix valued function P Xv that satisfies 
the following RH problem. 

RH problem for P Xv : 

(a) P Xv (z) is defined and analytic for z 6 Ug 0iXv \ £ for some <5o > <5- 

(b) -Pc^(z) satisfies the following jump relations on S° n Us >Xv ' 

p x ^ + {z) = p x ^{z)^ w{z) _^ {z) _ 2n °X forzG(S°nC ± )nC/ 5 ^, (4.1) 

P Xv>+ (x) = P Xv ,-(x) ( * x) _ x W{ ^y forxG(S°n(-l,l))n^. (4.2) 

(c) On dU$ jXv we have, as n ->oo 

P x Jz)N- 1 (z) =1 + 0(l/n), uniformly for z G 0Z7f lX „ \ S. (4.3) 

(d) For Ajy < 0, P Xv {z) has the following behavior 
iU*) = C>f! as*-*,. (4.4) 



s i |z - x v \ 2Xv J ' 

For Aj/ > 0, P Xl/ (z) has the following behavior 



O ( , „ ) , as z — ► x u from outside the lens, 



i r 
i iy 

1 2: I — -j^ 

O [ 1 i-2A y ) ' as 2 ^ from inside the lens. 

■2 X»y I 1 t 



(4.5) 



We will work as follows. First, we construct a matrix valued function that satisfies conditions 
(a), (b) and (d) of the RH problem for P Xv . For this purpose, we will transform (in Section [Q| ) 
this RH problem into a RH problem for P X J with constant jump matrices and construct (in 



Section 4.2) a solution of the RH problem for Pj„ . Afterwards, we will also consider (in Section 



4.3) the matching condition (c) of the RH problem for P Xv . 



4.1 Transformation to a RH problem with constant jump matrices 

Since h is analytic in U with positive real part, the scalar function 

v— 1 no 

W Xu {z) = (l-z) a / 2 (l + z)^ 2 h 1 / 2 (z)l[(z-x k )^ TJ 

k=l l=u+l 



(z-x u ) x », for z E (K Xi/ C\U)\ 
for z G (K r Xu n U) \ 



[Xu - z 



\ Ay 



(4.6) 



is defined and analytic for z € U \ (M U T Xv ). Here, we recall that K Xv and K Xv are the sets of 
all points on the left, respectively right, of T Xv . We seek P Xv in the form 

P,M = E n ^(z)P^(z)W x „(z)-°^(z)-^, (4.7) 
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where the matrix valued function E n>Xu is analytic in a neighborhood of Us, Xu i and E njXu will 



be determined (in Section 4.3) so that the matching condition (c) of the RH problem for P Xu is 
satisfied. 

Since P Xv has jumps on S n U^ Xu , and since W Xv has a jump on T Xu n U, the matrix valued 
function P x „ has jumps on the contour 

e Xv = (sur^jnt/^, 

see Figure |3|. The contour H Xu consists of 8 parts, which we denote by Si, ... , Eg, as shown in 
Figure [|. We write 

S° = E Xv \ {x v }, and T,° k = T, k \{x v }, for k = 1, . . . , 8. 




Figure 4: The contour S^. Here, S3 U S7 is the part of T Xi/ in Ug Xl/ , and the remainder is the 
part of S in U$^ Xv . 

In order to determine the jump matrices for P X J we need some information about the scalar 
function W Xv . Write 

Kl = Kl nc+, Kl l =Kl nc+, 

III _ isl rMT ZVIV 



K' u = KL. n 



K' w =K r x n ■ 



So, the sets , • • • , divide the complex plane into four regions divided by the real axis 
and the contour T Xi/ . By ( fTI)) and (|4~6|) , we have for 2 G fl K l x ) n f7, 



w{z)e 



-2*i\ % if zeKl u UK™, 
if z G if" U ifl v . 



Here, we recall that AZ n = i^i , , 

7 X J-UQ + l 



w(z)e 27riA % 

^. From this we see that 
W Xv) +(x)W Xv -(x) = w(x), for x G (x v -i,x u+ i) \ {x u }. 
By (^) we have on the contour T Xv , 

W Xv) +{z)W Xv) -{z)- 1 = e x ^\ for zeT^H U. 



(4.8) 



(4.c 



(4.10) 
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We now have enough information about W Xv to determine the jump matrices for PxJ . First, 



we determine the jump matrix on the lips of the lens. By (|4.1|), ( |4.7| ) and (4.8) the matrix valued 
function pj$ should satisfy on S° U S° U S° U the jump relation 

= P *IU Z ) ( w ( z )-i W l(z) l) 

= P«_( z ) ( e±2 ^ (4.11) 

where in e =t27nAl ' the + sign holds for z G £4 U Eg and the — sign for z 6 SjUSj!. 

Next, we determine the jump matrix on the interval. For x G S° U Eg we have by fl4.2|) , 
(O), (O), and the fact that ip + (x)ip^{x) = 1, 



= ^!-(^)^,-(^)- <T ^-(x)- nCT3 (_ w( ° ) _ 1 w{ x) ^w Xv , + (xr^ + (xr^ 



P (i) ° 



,d) ^ o ^(^w-^.+^r 1 ^,-^)- 1 



And finally, we determine the jump matrix on the contour that goes vertically through x v . 



For z G Eg U we have by (gj) and (|410|) 



= i'tWV^^+W' 4 = ^-O^™ 3 - (4.13) 

We then see that we must look for a matrix valued function P^} that satisfies the following 
RH problem. 

RH problem for Pi 1 ); 

(a) pJ^(z) is defined and analytic for z G Us 0jXv \ (E U r^) for some 5o > 

(b) Px\) (z) satisfies the following jump relations on S°^: 

P® + (x) = P£},-(x) jY forxGS?UEg, (4.14) 

pW + (z) = P^-(z)( e J nau j), far*€E§UEg, (4.15) 
pW + (z) = P^_(z)e x »™3, forz<EE§UE?, (4.16) 
PS,M = P^-( z )(jiX» J)' forzGS°US°. (4.17) 
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Figure 5: The contour T,q,. 



(c) For X v < 0, Pxl^ (z) has the following behavior 



P&H*) = 



z - x t 
z — X, 



| Z %ls 

Z Xj. 



clS Z > X]/ . 



(4.18) 



For X u > 0, pJ„ (z) has the following behavior 



o 



o 



\z 


Xy 


Xu 


\z - 


| 


-K 


\z 


Xy 


K 


\z - 


— *Ev | 


-K 


z - 


— Xy\ 


-K 


\z 


Xy 


-A, 


z - 


— Xy\ 


-K 


\z 




-A, 



cLS Z ^ X i/ from outside the lens, 



. cLS Z ^ X i/ from inside the lens. 



(4.19) 



Remark 4.1 Condition (c) follows from condition (d) of the RH problem for P x „, since 

P*A*) = E n ^(z)p£(z)W x „(z)-°^(zr n ^, 

where tp(z) is bounded and bounded away from near z = x u , and where W Xv {z) behaves like 
c\z — x v \^" as z — > x u , with a non-zero constant c. 

4.2 Construction of PxJ 

The construction of pj*' is based upon an auxiliary RH problem for in the C-plane with 
jumps on the contour = U^ =1 Ti consisting of eight straight rays, oriented as in Figure || 
which divides the complex plane into eight regions I- VIII, also shown in Figure |B[ We let 
A > -1/2. 

RH problem for *&\: 

(a) x I'a(C) is analytic for (£C \ S^. 
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(b) *a(C) satisfies the following jump relations on E^: 

*A,+ (C) = ^-(C)^ J), forCGTxUrs 



*A,+ (C) 
*A,+ (C) 
*A,+ (C) 



*xAo( e -Lx J). forcer 2 ur 6 , 
*A,-(C)e A ^ a , forCGr 3 ur 7 , 
*A,-(C)( e2 iiA J). forCGr 4 ur 8 . 



(4.20) 

(4.21) 
(4.22) 
(4.23) 



(c) For A < 0, ^a(C) has the following behavior as £ — ► 0: 

as C 0. 



•*(o-oi]<ji j<jr„ 



For A > 0, ^a(C) nas the following behavior as £ — ► 0: 



(4.24) 



*A„(C) = < 



of^x ^_aV ™c -O Jor^e [I. Ml. VI. \ ll. 



^( j^j-A j^j-A N ) ■ ™c -0 Tor c, £ I . IV. V. VIII. 



(4.25) 



We construct a solution ^a of this RH problem out of the modified Bessel functions I\±i/2 
and K^±i/2, and out of the Hankel functions #^2i/2 ano " ^A±i/2" ^ or C G I, we define VPa(C) d Y 



*a(0 = ^C 1/2 

For C G II, by 



A 7T A 77 . 



(4.26) 



*a(C) = 
For C G III, by 

*A(C) = 
For C G IV, by 



(4.27) 



-^C 1/2 / A -i(Ce- 



-^^^(Ce-- 



e |A7ri(T3 



(4.28) 



*a(o = Im-o 1/2 



<4(-o -<,(-c^ 



(A+j)™ 3 



(4.29) 
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For C G V, by 



*a(C) = o^(-C) 1/2 



iff 



(1) 



-(A+i)7Ti(T3 



-0 



(4.30) 




(4.31) 



-^C 1/2 / A+| (Ce^) - 
And finally, for £ G VIII, we define it by 



(4.32) 



1 



*a(0 = ^C 1/2 



/-iff 2,(0 -^?2i(c) 



A+ 



iff ±J !(c) -ffr ; i(0 



(2) 



3 (A+t)™ 3 



(4.33) 



Theorem 4.2 TTie matrix valued function ^f\, defined by ( 4.26| )- (|4.33| ), is a solution of the RH 
problem for 

Proof. The functions I\±x/2, ff\±i/2> -^l±l/2 an< ^ ^l±l/2 are defined and analytic in the com- 
plex plane with a branch cut along the negative real axis. So, the matrix valued function ty\ 
defined by ( 4.26|) - ([4.33|) is analytic in the respective regions, and condition (a) of the RH prob- 
lem is therefore satisfied. Condition (c) follows easily from |l], formulas 9.1.9, 9.6.7 and 9.6.9]. 
So, it remains to prove that jump conditions flPcD-flP3|) are satisfied. 

Jump conditions ( |4.20| ) and (4.22): By inspection, it is easy to see that these jump conditions 
are satisfied. 

Jump condition ( |4.21 ) for £ G r 2 : We use (4.27) to evaluate v I / a,+ (C) anci ( [4.26j ) to evaluate 
^A,-(C)- From ( [4.26D and [||, formulas 9.1.3, 9.1.4 and 9.6.3], the 1,1-entry and the 2,1-entry on 
the right of ( 4.21 ) are equal to 

*A,ii,-(C) + e- 2 ^ A , 12 ,_(C) 



2 2 A +2 

v^C 1/2 e- (A+ ^ m J A+i (C) 



v^C 1/2 / A+ l(Ce-^)e^^, 



(4.34) 



and 



*A,21,-(C) + e- 2 " iA *A,22,-(C) 



A — 7 



^C 1/2 e- (A+ ^"V A _i(C) 



(4.35) 
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respectively. By (|4.27[) we then see that the first columns of both sides of ( (4.21|) agree. From 
( 4.26|) , ( |4.27| ) and [l], formula 9.6.4], the second columns of both sides of ( |4.21| ) agree as well. 

Jump condition (|4.21 ) for ( G T 6 : We use (|4.30 ) to evaluate ^a,+ (C) an d (^-31| ) to evaluate 
v J / a,-(0- Since —Q = (,e m , we have, from ( 4 . 3 1| ) and formula 9.6.4], that the 1,2-entry and 
the 2,2-entry on the right of (4.21) are equal to 



A, 12, 



(0 = -^C 1/2 e^ A+| (Ce-) = -^(-C) 1/2 H^(-C)e^ 



and 



A,22, 



(0 = -^C 1/2 e^ M (Ce^) = ^(-0 1/2 <\(-C)e( A+ ^-, 



(4.36) 



(4.37) 



respectively. So, by ( 4.3C| ) we see that the second columns of both sides of ( |4.21[) agree. Since 
-( = (e™ we have, from ( pT3ll ), ( pL36] ), ( pT37|) a nd [g formulas 9.1.3, 9.1.4 and 9.6.3], that the 
1,1-entry and the 2,1-entry on the right of (|4.2l|) are equal to 



and 



^A,n,-(C) + e-^ A ^A,i2,-(C) 
= -v^(-C) 1/2 e- (A+ >V A+i (-C) + \^(-O l/2 e- (X+ ^Hf 

^A,2i,-(C) + e- 2 ^ A * A , 22 ,_(C) 
= ^M-C) 1/2 e-^I x ^((e^) - ^(-C) 1/2 e- (A+ ^-<\(-C) 

= v^(-C) 1/2 e-( A+ ^-J A _i(-C) - ^(-C) 1/2 e-( A+ ^ m <\(-C) 

2 A ~2 



(4.38) 



(4.39) 



respectively. By ( 4.30| ) we then see that the first columns of both sides of (4.21) agree as well. 
We now have proven that jump condition ( 4.21| ) is satisfied. 

Jump condition ( |4.23| ) : Similarly, we can prove that this jump condition is also satisfied. 
Here, we also use |IJ formula 9.1.35], and the details are left to the reader. This implies that 
the theorem is proved. □ 

>(i) 



Now, we explain how we get PxJ out of the solution ty\ v (depending on the parameter X u ) 

(4.40) 



of the RH problem for 1 §>\ v . We make use of the following scalar function, 



i\ogip{z) - i\ogtp + {x u ), 
-i\ogip{z) - i\og(p + (x v ), 



for Imz > 0, 
for Imz < 0, 
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Figure 6: The conformal mapping f Xv . For every k = 1, ... 8, the contour is mapped onto 
the part of the corresponding ray in f Xl/ (Us ;Xv ). 

which is defined and analytic for z G C \ M. For x G (— 1, 1) we have, since ip+(x)(f-(x) = 1, 
that f Xu ^-\-(x) = f Xv -(x), so that f Xv is also analytic across the interval (—1,1). The behavior 
of f Xu near x v is 



f Xv {z) = -j=L=(z -x u )+0 {{z - x v f) , 



s/l-.r 



, , as z 

2 
v 



Since f Xu is analytic near x v , and since f'(x u ) ^ 0, the scalar function f Xv is a one-to-one 
conformal mapping on a neighborhood of x u . So, if we choose 5 > sufficiently small, f Xu is a 
one-to-one conformal mapping on Ug )Xu and the image of Ug jXv under the mapping ( = f Xu is 
convex. 

For x £ (—1, 1) we have by (4.40) that f Xv (x) = arccos(x y ) — arccosx. So, f Xv (x) is real for 



x G (—1,1). If x > x v we have f x „{x) > 0, and if x < x v we have f Xl/ (x) < 0. Since f Xu is a 
conformal mapping, this implies that f Xu maps U$ :Xi/ nC + one-to-one onto f Xu (JJs jXv ) nC + , and 
Us, Xv n C_ one-to-one onto f Xv (Us jXv ) n C_. 

We now come back to the special choice of the contour T Xv , which we used to continuate 
our weight analytically, see Section |3.2| . For z G T Xv n C + we have aig (p(z) = arccosx^, by 
construction of T Xl/ , and for z G T Xv D C_ we have arg(/j(z) = — arccosav By ( 4.4C| ) we then 



have Ref x „(z) = 0, for z G r xv . This implies that the image of the contour Y Xv under the 
mapping £ = f Xv is the imaginary axis, which explains our choice of T Xv . 

We remember that the contour T, Xv was not yet completely defined. Now, we define the 
contours X^U^USgUEs as ^ ne preimages of the parts of the corresponding rays r2Ur4Ur6UF8 
in f Xv (Us >Xv ) under the mapping £ = f Xv (z), see Figure ||. We then have immediately that we 
can define 

P x lH*) = *\AnfxAz)), (4-41) 
and P x \) will solve the RH problem for P$ . 

Remark 4.3 We can use any one-to-one conformal mapping on U$ Xv to construct P x ] out of 



^x„- However, we have to choose it so as to compensate for the factor (f(z) n<J:i in (|4.7| ), In the 
next section we will see that our choice of conformal mapping will do the job. 

4.3 Construction of E n>Xu 

We recall that for every matrix valued function E nXv analytic in a neighborhood of Ug iXu , the 
matrix valued function P Xv given by 

P* V (A = E^ Xv {z)^ K {nf Xv {z))W Xv {z)~^^{z)~ n ^, (4.42) 
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satisfies conditions (a), (b) and (d) of the RH problem for P Xu . In this section we want to 
determine E n<Xv so that the matching condition (c) is satisfied as well. To this end we need 
to know the asymptotic behavior of at infinity, and use this to determine E n ^ Xv . At the 
end of this section we also show that E nXv is analytic in a neighborhood of U& Xu , so that the 
parametrix P Xv is completely defined. 

In order to determine the asymptotic behavior of at infinity, we insert the behavior of 
the Bessel functions at infinity into the matrix valued function ^i\ v , given by ( |4.26j )— (4.35). See 
Q formulas 9.7.1-9.7.4] for the behavior of the modified Bessel functions at infinity, and [|, for- 
mulas 9.2.7-9.2.10] for the behavior of the Hankel functions at infinity. Then, a straightforward 
calculation gives us the asymptotic behavior of at infinity. The behavior is different in each 
quadrant. For the upper half-plane we find as £ — > oo, 



**"<<> = 71 (-« T 



I + O 



uniformly for ( in the first quadrant, and 



/ + o ( i 



e f ff 3g-«(l73 e -|A„lW3 



uniformly for £ in the second quadrant. For the lower half-plane we find as £ — > oo, 



*A„(C) 



1 



V2 \~i 1 



I + O 



uniformly for ( in the third quadrant, and 



<M0 



i + o [ J 



1 1 



1 1 



(4.43) 



(4.44) 



(4.45) 



(4.46) 



y/2 \~i 1 

uniformly for £ in the fourth quadrant. 

Now, we use the asymptotic behavior ( [4.43| )-( |4,46| ) of ty\ v at infinity to determine E nfXv . 
We explain this only for the region 8U SjXv nKl v . The other cases are similar and the details are 
left to the reader. For z G dU$ >Xu n K x we have, since f Xu is a one-to-one conformal mapping 
on Ug jXv i that nf Xv (z) lies in the first quadrant, cf. Figure ||. So, we may use ( [4.43) to evaluate 
the asymptotic behavior of ^x u [nf Xv (z)) as n — > oo. Since Imz > 0, we have by (14.40 ), 



Using fljUp and ( gg ) we then find 



I + O 



■n 



x eX CT 3 (/ j + ( Xy )-^3 e -iA,7ri CT 3 T y^( z )-^ A r-i(- z ) ) 

as n — > oo, uniformly for z G dU$^ Xv n . So, in order that the matching condition is satisfied 
we define for 2 G ?7 D i£; , 



E n , x „(z) = N{z)W Xv {zY*e^™^+{xv) nai e-^ 



J_ (\ i 

V2\i 1 



(4.47) 
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Remark 4.4 With this E n ^ Xv we see that 



P x „(z)N-\z) = N{z)W Xv {zY*y + {x v ) 



11(73 



i + o\ - 

n 



x V+ (x u )- n ^W x „{z)-^N- L {z), 

as n — > oo, uniformly for z G 9Us, Xl , n ^ . Since l^+^j,)! = 1, and since W Xv as well as all 
entries of N are bounded and bounded away from on dUs sXl/ , the matching condition is satisfied 
on dU s <x v nKl v . 

Similarly, we use ( 4.40 ) and ( [4.42| ) together with the asymptotic behavior fl4.44|) -( |4.46| ) of 

at infinity to determine E n ^ Xu in the other regions. Straightforward calculations then show that 

we have to define E n ^ Xv {z) for z G U \ (M U T Xv ) as 



E ntXv {z) = E v (z)v + (x u r^e-^-L Q \ 



(4.48) 



where the matrix valued function E v {z) does not depend on n, is analytic for z £ XJ \ (R U 
and given by 



A^z)!^^)* 73 ^"™ 3 , for zEUHK^, 
N{z)W Xv (z) as e-i Kma3 , for z G [/ n 



N{z)W Xv {z] 



<?3 



1 

-1 



e 2 



iAj/7ricr3 



£„(z) = N{z)W Xv {z) c 







i or 



for z g £/ n if™, 



for z G f/n^ v . 



(4.49) 
(4.50) 

(4.51) 
(4.52) 



Now, everything is fine, except for the fact that E n>Xv is analytic in U \ (R U T Xv ), but we 
want it to be analytic in a full neighborhood of x u . This will be proven in the next proposition. 



Proposition 4.5 The matrix valued function E n>Xu defined by ( 4.48| )-( [4.52| ) is analytic in U \ 
((-oo,av_i] U [x^+1,00)). 



Proof. By (|148D it suffices to prove that E u is analytic in U \ ((— oo, x u ^\\ U oo)). We 

will check that E v has no jumps on {x v ^\,x v+ i) \ {x u } and (F Xi/ nU)\ {x„}, and in addition 
that the isolated singularity of E v at x v is removable. Let (x v -\, x v +\) be oriented from the left 
to the right, and let Y Xv n C/ be oriented so that it points away from x u , cf. Figure |3[ 

For x G (x u _i,x u ) we use (|4.50 ) to evaluate E v ^ + {x) and ([4.51| ) to evaluate E v _(x). From 
||[we have W Xvi+ (x) = w{x) l / 2 e Kni and W Xu -{x) = w(x)^ 2 e ;w \ Therefore, by (|3T5| ), 



(14501) and fl45p , 

^, + (x) = iV_(x) 



-w{x)~ 



w{x) 



w{x) 



CT 3 / 2 p A„ 7rio- 3 — i A„7ri.o- 3 



^,-(x). 



Hence, -Ey is analytic across (x v -±,x u ). Similarly, we have by (p^|) , (fO|, (jO|) and 
that is analytic across (x^rc^+i) as well. 
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For z £ (T Xv n U) n C+ we use (|4.50| ) to evaluate £^+(2) and Q4.49D to evaluate £v,_(,z 



From ( gig ) we have W^+^W^ = e*"**. Therefore, by ( fL49| ) and ( |450| ) 

= E V) - (z)e~^ m ^ W Xv - {z)~^ W Xu>+ {zY* e~ ^ mr7s 
= E v ,-(z), 



so that E v is analytic across (T Xu fl J7) Pi C + . Similarly, we have by ( 4.10[ ), ( f4.51| ) and ( 4.52| ) 



that E v is analytic across (T XiJ C\U)f\ C_ as well. We thus have proven that E v is analytic in 
U \ ((-oo, ^_a] U [x^+i, oo) U {x u }). 

It remains to prove that the isolated singularity of E v at x v is removable. We have by ( 3.17| ) 
and fl4.6|) that -D(z) behaves like c\\z — x u \ Xv and W Xv {z) like C2 1 z — x v \ v as z —* i v , where ci 
and C2 are non-zero constants. Therefore 



So, by ( |3.18[ ), each entry of N{z)W Xv {z) a ' i remains bounded as z — > av This implies by ( 4.49 ) 



(4.52) that each entry of E^, remains bounded J IS £ ^ Xjy, SO mat the isolated singularity of E v 



at Xj, is removable. Therefore, the proposition is proved. □ 
This ends the construction of the local parametrix near x v . 



We recall that we also wanted the local parametrix P Xv to be invertible, see Remark |3.5| . We 
will show that 

detP Xv = l. (4.53) 

This is analogous as in Section 7] and we will just give a sketch of the proof. Since E n>Xv is a 
product of four matrices all with determinant 1, see fl4,48|) -( fi.52 ), it suffices to prove from ( 4.42j ) 



that det^^ = 1. Using part (b) of the RH problem for we find that det is analytic 
in C \ {0}. If we then use the behavior of near stated in part (c) of the RH problem the 
isolated singularity of det ^l\ v at has to be removable, so that det ^!\ v is an entire function. 
Using the asymptotics of at infinity given by ( 4.43j )- fl4.46| ) we have that det 1 I / a„(C) — >• 1 as 



Q — > oo. By Liouville's theorem we then have that det ^!\ v = 1, so that also det P Xv = 1. 

5 Asymptotics of the recurrence coefficients 

In this section we will determine a complete asymptotic expansion of the recurrence coefficients 
a n and b n as n — * oo. Recall that a n and b n have been formulated in terms of the solution of 
the RH problem for Y, see ( |2.7| ) and ( |2.8j ). The asymptotic analysis of the RH problem for Y 
has been done in Section 3, and unfolding the series of transformations Y t— ► T \— ► S i— ► R, see 



[16, Section 9] for details, we find 



a n = J™. ( ~~2i~ + zR i2{z;n,w)J ^zR 2 i(z;n,w) + ^ 2 ), (5.1) 



->OC 



and 



b n = lim (z — zRu(z;n + l,w)R22(z;n,w)) . (5.2) 



->oc 
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Remark 5.1 We note that, see Lemma 8.3], 



R{z) - 1\\ = O -p-, , as n ^ oo, (5.3) 



n\z\ 

uniformly for \z\ > 2, where || • || is any matrix norm. Inserting this into ( |5.1| ) and ( |5.2[ ) we find 
the known asymptotic behavior of the recurrence coefficients, cf. [jDJ, 

a n = - + 0(l/n), b n = 0(l/n), as n — > oo. 

In the rest of the paper we will develop the 0(l/n) terms into complete asymptotic expansions 
in powers of 1/n. 

In order to determine a complete asymptotic expansion of a n and b n we will work as follows. 



In Section 5.1 , we will determine a complete asymptotic expansion of the jump matrix for R in 
powers of 1/n as n — » oo. As a result, we obtain in Section |5,2| a complete asymptotic expansion 
of R. The coefficients in this expansion can be calculated explicitly via residue calculus, and we 



will determine the order 1/n term. Finally, in Section |5.3| we will use this to prove Theorem 1.1 
5.1 Asymptotic expansion of A 

Denote the jump matrix for R as I + A. Then, from condition (c) of the RH problem for R, 
A(z) = P x „(z)N-\z)-I, lor zedU s , x „ and ^ = 0, . . . , n + 1, (5.4) 

A(z) = N(z) (^-i*^-*, °^N-Hz)-I, forze^Uu^+^J. (5.5) 
In this section we will show that A has an asymptotic expansion in powers of 1/n of the form 

a / \ A k (z,n) 

A ( z ) ~ Z_> k ' asrwoo, (5.6) 

k=l 

uniformly for z E and we will also determine the coefficients A^(z, n) explicitly. 

Remark 5.2 The n-dependance of the coefficients in the expansion will come from the factor 
(pj r {x v ) na ' i in the parametrices near the algebraic singularities x v . 

On the lips of the lens, A vanishes at an exponential rate, cf. |l6| , Section 7]. This implies 
for every k, 

A k (z,n) = 0, for z £ \ (U^c^J. (5.7) 

On the circles near ±1, the asymptotic expansion Q5.6| ) of A is known, see ]i~6| , Section 8]. The 
restriction of Afc(-,n) to dUs-i U dUs,i is given by |1(| (8.5) and (8.6)] and does not depend on 
n. It has a meromorphic continuation to Ug -i U Ug 0t i for some do > 5, with poles of order at 
most [(k + l)/2] at ±1. For details we refer to (l6|, Section 8]. 

So, it remains to determine the asymptotic expansion of A on the circles near the algebraic 
singularities. Fix u e {!,... ,n }. By (|4.42j ), ( [4.48|) and fy|), we have for z £ dUg tXv , 



A(z) = ^ + W"»e>i(! I 



x * Xu {nf Xv {z))y{zy™*W Xu {zy^N-\z) - I. (5. 
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Here, the matrix valued function *I'; Vj/ is constructed out of Bessel functions, which have a 
complete asymptotic expansion at infinity. This implies that ^ \ v (nf Xv (z)) also has a complete 
asymptotic expansion as n — > oo. Inserting the asymptotic expansions of the modified Bessel 
functions at infinity ||, formulas 9.7.1-9.7.4] into ( [4.26| ), and the asymptotic expansions of the 
Hankel functions [jl], formulas 9.2.7-9.2.10] into ( |4.27| ), we obtain 

( (-l) k s\ v ,k ~it\ v ,k 



*A„(C) 



1 



V2 \-i 1 



1 ok+i 



fc=i 



(5- 



as C oo, uniformly for £ in the first quadrant. Here, the constants s\ V) t and t\ v k are given 

by 



s\„,k = (K + ^,k) + (A„ - i, fe), = (A,/ + i, k) - (A„ - ^, k), 



(5.10) 



where 



(«/, fc) 



(4i/ 2 - 1)(4^ 2 - 9) . . . (4z^ 2 - (2k - l) 2 ) 
2^! ' 



For z S dU$ jXv fl if* we have, since / xv is a one-to-one conformal mapping on U$ jXv that nf Xu (z) 
lies in the first quadrant, see Figure ||. So, we may use (|5.9| ) to determine the asymptotic 
expansion of ^Ai/fa/avC 2 )) as in oo. Since Imz > we have by ( [4.40 ), 

e -in/.„(*) = ¥ , + ( Sl/ )-n v ,( z )«. 



Therefore, by (5.9), 

_L_ fl i 

V2 



n<T3 



^U^zfn* \i{-l)H Xv>k s Xv , k 



e J 



<T3 



as n — ► oo, uniformly for z G dU$ tXv D if^- Inserting this into (5.8) we have by ( f4,49j ), and the 
fact that v? + (x l ,) n remains bounded and bounded away from as n — > oo (since |<£ + (xj,)| = 1), 



A(*) 



: E u (z)<p+ (x u ) 



-l) fe SA„, fc 



-(-l) fc *A„,fc S\„,k 



(5.11) 



as n — > oo, uniformly for z 6 dU$ >Xu D iCj . Similarly, we find the same asymptotic expansion 
on the other regions of dUg Xl/ . The details are left to the reader. Thus, for z £ 9Ug Xi/ the 
coefficients of the expansion ([)!]) for A are given by 



Ak(z, n) 



2 k+1 f Xv (z) k 



\U(T3 



-(-l) k t\ v ,k S Xv ,k , 



V+{x v )- n ^E-\z). (5.12) 



Remark 5.3 These coefficients depend on n through the factors (p+(x u ) zizncr3 . Since |(^+(xj,)| = 
1, the coefficients A k (z,n) remain bounded asm oo, which is necessary to get an asymptotic 
expansion of the form ( |5.6| ). 
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We note that /~ is analytic in C \ ((— oo, —1] U [1, oo)) except for a pole of order k at x v , 
see the discussion at the end of Section [4^. Prom the proof of Proposition 4.5 and the fact that 
det E u = 1 we have that E v as well as E~ x are analytic in U \ ((— oo, :r„-i] U °°))- So, the 

restriction of Afc(-,n) to dUs jX „ has a meromorphic continuation to a neighborhood Us 0:Xl , of z v 
for some 5q > 5, with a pole of order k at x u . 



5.2 Asymptotic expansion of R 



We recall that A possesses an asymptotic expansion in powers of 1/n of the form ( |5.6p with 
oscillatory terms in the expansion. Following the argument that leads to jj], (4.115)], this implies 
that R itself possesses an asymptotic expansion in powers of 1/n given by 



R(z; n,w) ~ I + 



Rk(z,n) 



as n 



k=l 



oo, 



uniformly for z £ C \ Ejj, Here, for every k and n, 
i? fc (-,n) is analytic in C \ (U^dC/^J, 

and 

tfitO^n) = 0(1/2), 



as z 



00. 



(5.13) 



(5.14) 



(5.15) 



The n-dependance in the coefficients Rk(z, n) arises through the oscillatory terms in the expan- 
sion of A. 

We will now determine, similar as in [[l(], Section 8], the coefficient Ri(z,n) explicitly. Ex- 
panding the jump relation R + = i?_(J + A), and collecting the terms with 1/n we have 



i?i,+(s,n) - Ri,-(s,n) = Ai(s,n), for s G U^dUg^ 



(5.16) 



which is, together with ( p. 14 ) and ( |5.15| ) an additive RH problem. This can easily be solved 
using the Sokhotskii-Plemelj formula, but in our case we can write down an explicit solution 
as follows. Since Ai(z, n) is analytic in neighborhoods of z = ±1 and z = x u for v = 1, . . . no, 
except for simple poles at those points, see Section 5.1, we can write 



and 



Ai(z,n) 



Ai(z,n) 



+ 0(1), as z->l, Ai(z,n) 



z+1 



+ 0(1), 



as z 



1, 



+ 0(1), 



as z 



n 



for certain constant matrices A^(n),B^(n) and Oi ( 

Remark 5.4 Since Ai(s,n) is independent of n for s E dUs-\ U dUg^i, see Section [Tl|, the 
residues ^4^^(n) and B^\n) of Ai(z,n) at z = 1 and z = —1, respectively, are also independent 
of n. The n-dependance of the residue C^\n) at z = x u follows from the oscillatory terms 
ip + (x u ) ±nr73 in Ai(s,n) near x u , see ( 5.12 ). 



By inspection we then see that 



AW {n) BW(n) 



Ri(z,n) 



z-1 

#(n) 
z-1 



+ 



z + 



"" C«(n) 



l/=l 

no 



1 ^ z 



Ci l \n) 



z + 



for 2 eC\(uS 1 f/^). 



Ai(z,n), fczeU^V. 



(5.17) 
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satisfies the additive RH problem ( 5.14| )- (|5.16| ). So, we need to determine the constant matrices 
A^\n),B^\n) and C { v\n) for v = 1, . . . ,uq. For AW(n) and B^\n) we have found, see p6| 
Section 8], 



4a 2 



16 



-D a ' A 

oo 



(5.18) 



16 00 



(5.19) 



which is clearly independent of n. It remains to determine the residue C^\n) of Ai(z,n) at 
2 = x u , for every i/ = 1, . . . , n . Fix v 6 {1, . . . , no}. Since E v (z) and E~ l (z) are analytic in a 
neighborhood of z = x„, and since 



fxA z ) 1 = \A - HO(l), as z ^ x,,, 

we have by ( p. 10 ) and (5.12) that the residue (n) of A\(z,n) at = x v is given by 



-2\i -2\ v 
2\ v 2\l 



?\v + {x u )-^E-\x u ). 



(5.20) 



We want to simplify this expression. So, we need to find convenient expressions for E v {x y ) 
and E~ l {x u ), and substitute these into ( 5.20 ). Since E u is analytic near x v we determine E v (x u ) 
by the following limit 

E v {x v ) = lim E v (x) = lim E Ut+ (x). 

Here, we take the limit from x to x v on the real axis from the right. The last equality follows 
from the fact that E v has no jumps on (x u , x u +i), see Proposition 45. From (|3.18| ) and Q4.49Q 
we then find 



E v {x v ) = lim D£ 

xlx v ;x<=:K 



aj r {x)+aj r {x) 1 aj r {x)—aj r {x) 1 

2 2l 

aj r (x)—a,j r (x)~ 1 a+(x)+a+(x) _1 

=2i 2 



DUx) 



(5.21) 



The Szego function satisfies D + (x) = y?w{x)e '^"W for x £ (x u ,x u+ i), see Lemma |3.4| , where 
ip u is given by (fjg) . By (fD| we have Wa^+fc) = vM^e - *"™, so that by Q and (|^) 



li m ^W+WglA^Tri = e #^(a;„) e -iA„7rj = g-^*^ 



Inserting this into ( |5.21 ) and using the following identities, which hold for x £ (—1, 1), 

a+(x) + a + (x)- 1 _ e~x /2 
2 ~ ^(l-x 2 ) 1 / 4 ' 

a+(x) - a+(x) _1 _ e'x . _ ^_i/ 2 



2?: 



v^l-x 2 ) 1 ^ 



we then find 

E v (x v ) 



e 4 



^(1-^)1/4 



(^)V2 i^ + (x„)-V2 



V?+(X^ 



,1/2 



e 2 



(5.22) 
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Taking inverse, we find 



F -l( r )- 6 1 A* V «TB ( y+(^) 1/2 -^+(^) V2 \ D -,T, , V> ,,, 



Now, we insert Q5.22 ) and ( 5.23j ) into ( |5.20| ). Using the identity tp+(x v ) = exp(i arccos x v ) 
we then find after a straightforward calculation that the residue C{/ '(n) of Ai(z, ra) at z = Xj, 
is given by 

C «W = ^(^;<:) ^W', (5-24) 



where 



C^,n( n ) = ~2^ Xv + Sin arccos Xu ~ (J>1/ )' (5.25) 
C u ,i2(n) = % -)? v - % -\ v Xv sin(2n arccos x v - 



2 " 2 

~ 2 



Aj, \J\ — xl cos (2n arccos x v — $^ ) , (5.26) 



Cv,2i(n) = ^A^ - ^AyX I/ sin(2n arccos - <Jv) 



+ -A^\/r^l^cos(2n arccos a;,, — <&,,). (5.27) 

This ends the determination of Ri(z,n). 

For general fe, we get that Rk(z,n) in the region C\ (U^J t/j^J is a rational function with 
poles at ±1 and at the algebraic singularities x u . The residues at 1 and — 1 are denoted by 
A^ k \n) and B^ k \n) respectively and may depend on n. The residue at every x v depends on n 
and is denoted by cj) k \n). We then get 

„, . ,4( fc )( n ) B (k)(„) J^ fc) (n) „„ /2 , 

R k (z,n) = Y + -V + V i; +Q(1A 2 ), asz^oo. 

2 — 1 2+1 ^— f 

The residues A^{n), B^(n) and (n) can be determined in a similar fashion, but for our 
purpose it suffices to know R\{z, n). 

5.3 Proof of Theorem |T| 

We are now ready to determine a complete asymptotic expansion of the recurrence coefficients 
a n and b n . The idea is to insert the asymptotic expansion ( [5.13 ) of R into ([O]) and ( [5.2|) . 

Proof of Theorem p..l| . We recall that, see (5.1), 

a^= lim ( — -2^ + zR\2(z; n, w) J ( zRai{z; n, w) + 



2t "~"\~'~>-tJ y~H«xv~ ) ",-V 2iD 2 

We may take the limit z — > oo in the asymptotic expansion ( |5.13| ) of i?, cf. Section 9], to 
obtain 



, Dl , ^ A ^2 (n) + S w (n) + c w 2 (n) 



2i ^— ' n 

fc=i 



^ 4i>) + *ff(") + E£iO") i \ 
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as n — > oo. Expanding this we find a complete asymptotic expansion of o^, and this leads to 
a complete asymptotic expansion of a n in powers of 1/n as n — > oo. By ( |5.18|) , ( |5.19| ), ( |5.24| ), 
( 5.26 ), ( 5.27 ) and ( |5,28 ) the first terms in the asymptotic expansion of a 2 n are 



1 1 

4 + 2i 



D~J A${n) + B$(n) + £ 



,(1) 



!/=l 



,(1) 

!/,21 



/J 



l/=l 



7 - "tt "\A ~ ^ cos ( 2n arccos x v - — + O ( — ^ ] , 
4 ^— ' 2 v n \n z I 

u=l x 7 



(5.29) 



as n — * oo. From this we then get, after a simple calculation, that the coefficient with the 
1/n term in the asymptotic expansion of a n is given by ( |1.4| ). So, the statements about the 
recurrence coefficient a n are proved. 

Similarly, we can prove the statements about the recurrence coefficient b n . If we take in ( |5.2| ) 
the limit z — > oo in the asymptotic expansion (|l3|) of R, cf. H, Section 9], we find 



6 n ~ lim -z y~] 



( J R fc )ii(z,n + 1) ^(R k ) 2 2(z,n) 



\k=l 



(n + 1) 



l) fc + ^ 



fc=i 



rr 



fc=i 



4?(n + 1) + g^(n + 1) + yzu c;a(" + 1 

(n + l) fc 

AgJ(f»)+sgHn) + E^ciSi(n)\ 



+ 



n' x 



(5.30) 



as n -> oo. From this we get a complete asymptotic expansion of b n in powers of 1/n, and 
by fl538|) , (|5TT9D , QETgg ), ( ggg ) and flOGD the coefficient with the 1/n term in the asymptotic 
expansion of b n is given by 



Bi(n) 



A W („+!) + («))-( B^ (n + 1) + B $ (n) 



E( c Si( n + 1 )+ c . 



(i) 

i>,22 



i/=i 



sin^(2n + 1) arccos x u — Q u + arccos x u ^j 

— sin^(2n + 1) arccos x v — & u — arccos x l 
= — X u \/l — xf, cos((2n + 1) arccos x v — . 
Therefore, the theorem is proved. 



(5.31) 



□ 
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